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validity of this identity is discussed in detail. Most importantly, we show that due to the infrared 
finiteness of the gluon propagator, the QCD charge obtained with either definition freezes in the 
deep infrared, in agreement with theoretical and phenomenological expectations. 
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1. Introduction 

One of the most difficult problems in QCD is to understand the interface between the pertur- 
bative and non-perturbative regimes. Both sophisticated theoretical tools [nl |2|] as well as more 
phenomenologically oriented approaches [{3|] indicate that this connection is not abrupt, but rather 
smooth. 

Without any doubt, a continuous interpolation between the perturbative and the non-perturbative 
region is intimately related to the behavior of the QCD fundamental coupling [Q]. However, we 
know by now, that the smoothness in this transition can not be achieved with perturbative assump- 
tion of a coupling with singular growth (Landau pole) in the infrared. Indeed, all accumulated evi- 
dence points toward the need of a freezing of the QCD effective charge at small momenta ^ g ^jj, 
where a(q 2 ) develops an infrared fixed point and QCD has a conformal window at low energy [g]. 

The infrared finiteness of the effective charge can be considered as one of the manifestation 
of the phenomenon of dynamical gluon mass generation [ffl, |6|] revealing in this way, its profound 
connection with the most fundamental Green's functions of QCD, such as the gluon and ghost 
propagators [Q]. Indeed, the basic ingredients that enter in its definition must contain the right 
information and be combined in a very precise way in order to endow the effective charge with the 
required physical and field-theoretic properties [§]. 

In this talk, we will show how different QCD Green's functions can be combined in order to 
form renormalization group (RG) invariant quantities which eventually may be associated to a def- 
inition of an effective charge [Qj. Specifically, we will consider firstly the effective charge obtained 
within the pinch technique (PT) framework 0], and its correspondence [||] with the background- 
field method (BFM) The PT effective charge constitutes the most direct non-abelian gener- 
alization of the familiar concept of the QED effective charge. The second definition involves the 



ghost and gluon self-energies [10], in the Landau gauge, and in the kinematic configuration where 



the well-known Taylor non-renormalization theorem [11] becomes applicable. 

2. Definitions and ingredients 

Let us introduce some of the basic ingredients necessary for the definition of the two effective 
charges we want to study. In the Landau gauge, the full gluon propagator An V (q) is transverse, and 
omiting the color indices, its general form is given by 

A flv (q) = -iP^ v {q)A(q 2 ), with P„ v ( 9 ) = g „ v - ^ , (2.1) 

where the scalar function A(q 2 ) is related to the all order self-energy Un V (q) = P^ v (q)Yl(q 2 ) 
through A _1 (<7 2 ) = q +iTl(q 2 ). The full ghost propagator D(q 2 ) and its dressing function F(q ) 
are related by 

Dtf) = ^. (2.2) 
q 1 

In our construction, a special role is played by the auxiliary two-point function A^ v (q), repre- 
sented in Fig. [j] and defined as [12] 



A llv {q)=g llv G{q 2 ) + ^L{q 2 ) = -ig 2 C A jH$D(k + q)AP a (k)H av (k,q), (2.3) 
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where Ca is the Casimir eigenvalue of the adjoint representation [Ca = N for SU(N)], and 
J k = jJ. 2e (2n)^ d J d d k, with d = 4 — e the dimension of space-time. The vertex H^ v (k,q) is also 
represented in Fig. [j], and its tree-level counterpart is given Huty = ig^v An additional constraint 
on the behavior H^ v (k,q) is imposed by the WI (Ward identity) 

q v H llv (k,q) = -iT ll {k : q), (2.4) 

where T^(k,q) is the all-order ghost vertex, with k representing the momentum of the gluon and q 
the one of the anti-ghost; at tree-level T^\k,q) = — q^. 



H ai/ (k, q) = 




k + q 

Figure 1: Diagrammatic representation of H. 



2.1 The pinch technique effective charge 



The heart of the PT effective charge definition lies on the construction of a new effective gluon 
propagator, A(q 2 ), which captures the running of the QCD /3 function, exactly as happens with 
the vacuum polarization in the case of QED [||, ^, §]. Already at one-loop level, the PT gluon 
propagator displays the desired coefficient in front of the perturbative logarithm, namely 



1+- 



(2.5) 



where b = 11Ca/487T 2 is the first coefficient of the QCD /3 -function when the number of fermions 
rif = (quarkless QCD). 

In addition, to all order A~ l (q 2 ) is universal (i.e. process-independent), and therefore it does 
not depend on the details of the process where it is embedded as shown in Fig. |2[ 






(6) (c) 
Figure 2: The universal PT coupling. 



One important point, explained in detail in the literature, is the (all-order) correspondence 
between the PT and the Feynman gauge of the BFM [M 0J. In fact, one can generalize the PT 
construction ^ in such a way as to reach diagrammatically any value of the gauge fixing parameter 
of the BFM, and in particular the Landau gauge. In what follows we will implicitly assume the 
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aforementioned generalization of the PT, given that the main identity we will use to relate the two 
effective charges is valid only in the Landau gauge. 

Due to the Abelian Wis satisfied by the PT- BFM Green's functions, the renormalization con- 
stants of the gauge-coupling and of the PT gluon self-energy, defined as 

g(ll 2 ) =Z- V)*o; A(^,At 2 ) =Z A V)A (<7 2 ), (2.6) 

where the "0" subscript indicates bare quantities, satisfy the QED-like relation 

Z,=Z A 1/2 . (2.7) 

Thus, it follows immediately that the product 

d (q 2 ) =g 2 Q Ao(q 2 ) =sV)A(<? 2 ,At 2 ) =d(q 2 ), (2.8) 

retains the same form before and after renormalization, i.e., it forms a RG-invariant {fL -independent) 
quantity [[I], [|, pt]. For asymptotically large momenta one may extract from d(q 2 ) a dimensionless 
quantity by writing, 

d(q 2 ) = ^P-, (2.9) 
T 

where g 2 (q 2 ) is the RG-invariant effective charge of QCD; at one-loop (use Eq. ( |23| ) into (^8|)) 

iV) = i , r-fjWr^ = / ! : x , (2-10) 



QCD 

where Aqcd denotes an RG-invariant mass scale of a few hundred MeV. 

Being a direct consequence of the Wis satisfied by the PT Green's function, Eq. (2.8) may be 
employed either perturbatively or non-perturbatively, provided that one has information on the IR 
behavior of the PT-BFM gluon propagator A(q ). 

However, thanks to a general relation connecting the PT-BFM A(q 2 ) and the conventional 
gluon propagator A{q 2 ), all the non-perturbative information we have gathered about A(q 2 ) may 
also be used. Specifically, the aforementioned formal all-order relation states that Jl3| ] 

A(q 2 )=[l+G(q 2 )] 2 A(q 2 ), (2.11) 

where G(q 2 ) is the form factor of the g^ v component appearing on the definition of A^ v given in 
Eq.(^3|). Note that, due to its BRST origin, the above relation must be preserved after renormaliza- 
tion. Specifically, denoting by Za the renormalization constant relating the bare and renormalized 
functions, Aq v and A^ v , through 

g» v +A^(q,H 2 )=Z A (n 2 )[g» v +AZ v (q)}, (2.12) 



then from (2.7) and (2.11) follows the additional relation 



Z g l =Z l / 2 Z K , (2.13) 



which it will be useful in the following subsection. 



4 



Non-perturbative Green's functions and the QCD effective charge 



Arlene C. Aguilar 



At lowest order, it is straightforward to verify, that the role of the function the 1 + G(q 2 



obtained from Eq. (|2.3|), is to restore the /3 function coefficient in front of UV logarithm. Explicitly, 
we have at one-loop (in the Landau gauge) jl2| ] 

13 C A g 2 



l+G(q z ) 



,9C A g 2 (q 2 
1 + 448F ln '^ 



q 



1 + 



2 487T 2 



In It 



(2.14) 



Using Eq. (2.1 1) we therefore recover the Ar l (q 2 ) of Eq. (E3b, as we should. 

Then, non-perturbatively, one substitutes into Eq. ( 2.11 ) the 1 + G(q 2 ) and A{q 2 ) obtained from 



either the lattice or SD analysis, to obtain A{q 2 ). As explained above, the combination formed by 



d(q 2 



g 2 A(q 2 



i2 ' 



[l+G(q 2 )X 

is independent of the renormalization point fL i.e. a RG-invariant quantity. 



(2.15) 



2.2 Gluon-ghost vertex 

Another possibility for defining the QCD effective charge can be obtained starting from the 
various QCD vertices. The basic idea behind is to recognize the RG-invariant quantities we may 
form out of these vertices. The downside of this construction lies in the fact that it involves all the 
momentum scales present in the vertex in question, and further assumptions about their kinematic 
configuration need to be introduced, in order to express the charge as a function of a single variable. 
The ghost-gluon vertex has been particularly popular in this context, especially in conjunction with 
Taylor's non-renormalization theorem and the corresponding kinematics Jl0|]. 



Taylor 
kinematics 



k + q / \ Q 

z' g = z x z A ll ' 2 z- 1 



/ \ 9 



ryl _ 7 -l/2 7 _l 

^ a — 



Figure 3: The ghost-gluon vertex and the Taylor kinematics. 
For the case of the ghost-gluon vertex (see Fig. 0), the renormalization constants involved are 



A(q 2 ,^)=Z A \ix 2 )A (q 2 ), 
r v (k,q^ 2 )=Z^ 2 )T v (k,q), 



F(q 2 ,^)=Z;\^ 2 )Fo(q 
g =Z g ,(ll 2 )g'. 



(2.16) 



Z\Z A 1 ^ 2 Z C 1 , does not have to coincide with the Z g 



Notice that a priori Z g > defined as Z g , 
appearing in ( |2.6| ); however, as we will see in the next section, they do coincide by virtue of the 
basic identity we will derive there. 

Choosing the special Taylor's kinematic configuration, where the incoming ghost has a van- 
ishing momentum (i.e. — > —q^), one may impose the following additional condition one the 



renormalization constant Z\ (valid only in the Landau gauge), namely Z\ = Z g iZ A 
which follows immediately that 

7-1 _ 7 !/2 7 



1/2 



Z, = 1, from 
(2.17) 
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Thus, the product 

?(q 2 ) =g' 2 A(q 2 ;Ll 2 )F 2 (q 2 ;H 2 ) = g' 2 A (q 2 )F 2 (q 2 ), (2.18) 
forms a dimensionful jx -independent combination. Therefore, for asymptotically large q 2 , in anal- 



ogy to Eq. (2.9), one can define an alternative QCD running coupling as 



Using then Eq. ( |2.14| ), and the fact that 



F-\q 2 ^ 



1 + 



9C A g 2 



4 48?r 



In 



(2.19) 



(2.20) 



it is straightforward to verify that g 2 h (q 2 ) and g 2 (q 2 ) display the same one-loop behavior, since, 
perturbatively the function \ + G{q 2 ) is the inverse of the ghost dressing function F(q 2 ). As we will 
see in the next section, this is nothing more than the one-loop manifestation of the more general 
identity relating G(q 2 ) and F(q 2 ). 

3. A special relation between Green's functions 

In this section, we sketch the main steps needed to derive the central identity, valid only in 
the Landau gauge, relating the ghost dressing function with a particular combination of the form- 



factors G(q ) and L(q ) appearing in the tensorial decomposition of A^ v given in Eq. (2.3) 



> O * ) 1 = ( 

q 



q 




k + q 

Figure 4: The SDE for the ghost. 
First, consider the standard SD equation for the ghost propagator (Fig ^|), 
iD \q 2 )=q 2 + ig 2 C A J^A^ (k)T v (k, q)D(k + q) . 
Then, contract both sides of the defining equation (23) by the combination q^q v to get 
[G(q 2 )+L(q 2 )]q 2 =g 2 C A f q p AP a (k)q v H (JV (k,q)D(k + q). 

J k 



(3.1) 



(3.2) 



Using Eq. (2.4) and the transversality of the full gluon propagator, we can see that the rhs of 
Eq. (3.2) is precisely the integral appearing in the ghost SDE ( p.l[ ). Therefore 



[G(q 2 )+L(q 2 )]q 2 =iD- i (q 2 )-q 2 , 
or, in terms of the ghost dressing function F(q 2 ) [viz- Eq. (^2|)] 



(3.3) 



\+G{q 2 )+L{q 2 )=F-\q 2 ) 



(3.4) 
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The above relation, derived here from the SD equations of the theory, has been first obtained 
in in the framework of the Batalin-Vilkovisky quantization formalism. As was shown there, 
the relation is a direct consequence of the fundamental BRST symmetry. 

Now, let us construct the dynamical equations governing the behavior of the functions G(q 2 ) 
and L(q 2 ). The tensorial projection of both functions in terms of A^ v may be obtained from 
Eq. (^3|), where in d dimensions, we have 



G{q 2 ) 



1 



(^AjJ-gYA^v), L(q 2 



1 



(d-l)q 2 

which then gives, in terms of the SDE integrals 



(d-l)q 2 



(rf/^v-^), (3.5) 



G(q 2 ) 
L(q 2 ) 



g 2 c A 

d-l 

8 2 C A 
d-l 



AP a (k)H ap {k,q)D(k + q)+i\ f q p A p(J {k) T a (k,q)D{k + q) 

q 1 Jk 

i^jqPA pa (k) r ff (k, q)D(k + q) +^ A" ff (*) H ap (k, q)D(k + q) 



(3.6) 



In this point some words about the renormalization and approximations we will employ in 
Eq. Q6) are in order. Let us start with the renormalization procedure. As mentioned before, 
since the origin of ( 3A ) is the BRST symmetry, it should not be deformed after renormalization. 
Combining the definitions of ( 2.12 ) and ( 2.16D , we see that in order to preserve the relation (^4j) 
we must impose that Z A = Z c . In addition, by virtue of (2A), and for the same reason, we have that, 
in the Landau gauge, T v (k,q) and H av (k,q) must be renormalized by the same renormalization 
constant, namely Z\ [viz- Eq. ( |2.16 )]; for the Taylor kinematics, we have that Z\ = \. 

Then, approximating the two vertices, H llv (k,q) = ig^ v , and T^(k,q) = —q^, by their tree- 
level values, then, setting f(k,q) = (k- q) 2 /k 2 q 2 , one may show that [Q] 

F-\q 2 ) = Z c +g 2 C A f[l-f(k,q)}A(k)D(k + q), 

g 2 c A 



l+G(q 2 ) =Z C + 



d-l Jk 



[(d-2)+f(k,q)]A(k)D(k + q), 



L(q 2 



g 2 c A 

d-l.)k 



l-df(k,q)]A(k)D(k + q), 



(3.7) 



which clearly satisfies Eq. (p. 4]). 

We next go to the Euclidean space, by setting — q 2 = q\, and defining Ae(<7 e ) = — A(— g E ), 
De(^|) = — D(— q\), and for the integration measure J k = iJ kE . Then, suppressing the subscript 
"E" and setting q 2 = x, k 2 = y, we have that k-q = ^/xycos 6, and so (k ■ q) 2 jq 2 = jcos 2 6, and 
(k + q) 2 = x + y + 2^/xycosd, we arrive at (see details in [Q]) 



1 + G(jc) 



L(x) 



F'Hx) 



&sC A 

I6n 



F(x) 



WsC A 

12% 



Fix) 



/o^K 3+ s) Aw+ r^( 3+ i) AMFM 

A(y)F(y) 



dyy z A(y) +x dy 



MsC A 
167T 



Fix] 



J*dy y(3-l) A(y) + £dy (3 - -) A(y)F(y) 



(3.8) 
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Then, it is easy to see (e.g., by means of the change of variables y = zx) that if A and F are IR 
finite, Eq. ( |3.8| ) yields the important result L(0) = Let us now assume that the renormaliza- 
tion condition for F(x) was chosen to be F{ji ) = 1. This condition, when inserted into the third 



equation of ( |3.8| ), allows one to express Z c as 



1 + 



167T 



M 2 JO 



A(y) + / dy ( 3 



A(y)F(y) 



(3.9) 



and may be used to cast ( ]3.8| ) into a manifestly renormalized form. Note that if one choses 
F(fL 2 ) = 1, then one cannot choose simultaneously G(fL 2 ) = 0, because that would violate the 
identity of Eq. (3.4), given that L(fL 2 ) ^ 0. In fact, once F(fL 2 ) = 1 has been imposed, the value of 
G(fL 2 ) is completely determined from its own equation, i.e. the first equation in (3.8). 

In addition, in the MOM scheme A(q 2 ) and A(q 2 ) cannot be made equal at the renormalization 
point, since Eq. Q2.11| ) implies A^(jl 2 ) = \i 2 [l + G 2 (/i 2 )] 2 . 



Now, let us to return to the couplings, and discuss the implications of the identity given by 
Eq. (3.4). First of all, comparing Eq. (2.13) and Eq. (2.17), it is clear that g(jx) = g'(fl), by virtue 



of Za = Z c . Therefore, using Eq. ( |2.11| ), and the definitions given in Eqs. ( |2.8| ) and Q2.18D , one can 



obtain a relation between the two RG-invariant quantities, r(q 2 ) and d(q 2 ), namely 



?(q 2 ) = [\+G(q 2 )YF 2 ( q 2 )d(q 2 ). 
From this last equality, follows that a VT (q 2 ) and a g h(q 2 ) are related by 

a gh (q 2 ) = [\ + G(q 2 )] 2 F 2 (q 2 )a PT (q 2 ), 



(3.10) 



(3.11) 



After using Eq. (3.4), we have that 



a„{q 2 ) = a gh ( 



1 + 



L(q 2 



l+G(q 2 ) 



(3.12) 



Evidently, the two couplings can only coincide at two points: (i) at q 2 = 0, where, due to the 
fact that L(0) = , we have that a g h(0) = ot PT (0), and (ii) at q 2 = °°, given that in the deep UV 
L(q 2 ) approaches a constant. Note in fact that the two effective charges cannot coincide at the 
renormalization point jj,, where a g h(/i 2 ) = [1 — L(/x 2 )] 2 a PX (/x 2 ); this can be understood also in 
terms of the discussion following Eq. (|3.9|). 



4. Schwinger-Dyson input and numerical analysis 

Now we are in the position to compute the QCD effective charges defined above, using as 
input the non-perturbative solutions of SDE for the various Green's functions appearing in their 
definitions. More specifically, we will solve numerically a system of three coupled non-linear 
integral equations in the Landau gauge, containing A(q 2 ),F(q 2 ), and G(q 2 ) as unknown quantities. 

Once solutions for these three functions have been obtained, then L{q 2 ) is fully determined by 



its corresponding equation, namely the second one in Eq. (3.8). 
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Figure 5: The new SDE for the gluon one-loop dressed diagrams. 



The two SD equations determining F(q 2 ), G(q 2 ) are given in Eq. (3.8). The SD equation 
governing A(q 2 ), is given by [12] 



[1 +G(q 2 )] 2 A- l (q 2 )P^ v ( q ) = q% v (q) +/£» M v, 



(4.1) 



/=i 



where the diagrams (a;)^ are shown in Fig. ||| As explained in [12], due to the abelian WI sat- 
isfied by the fully-dressed vertices in the PT-BFM scheme, we have that q^[(ai)^ v + (02)^] = 
(f' [( a 3)/iv + {dA)tiv\ = 0- This last property enforces the the transversality of the gluon self-energy 
"order-by-order" in the dressed-loop expansion, which is one of the central features of the gauge- 
invariant SD truncation scheme defined within the PT-BFM framework. 

After introducing appropriate Ansatze for the aforementioned fully-dressed vertices, we finally 
arrive at the integral equation 



[1 + G(<7 2 )] 2 A-V) = I 1 



8 2 C A 



A(k)A(k + q)f y + A(k)f 2 



+ g 2 C A 





\k 2 {k ' q)1 ] 


.3 A 


q 2 _ 



k 2 {k + q) 2 
D(k)D(k + q)-2 J ' D(k) 



(4.2) 



with 



fx 
h 



20q 2 + \U 2 -6(k + q) 2 + 



(q 



2\2 



(k-qf 



27 
"2 



k 2 



+ 



q 



(k + q) 2 

(k-q) 2 



• + 



20 10 

¥ + q 2+ k 2 {k + q) 2 



+ 



2(k + q) 
q 2 k 2 



21 



(k + q) 2 (k + q) 2 k 2 (k + q) 2 q 2 (k + q) 2 



(4.3) 



The important point is that, by virtue of the poles introduced into the equation through the particular 
Ansatze employed [111 |TJ[ O], one obtains an IR finite solution for the gluon propagator, i.e.a 
solution with A -1 (0) > 0, in complete agreement with a large body of lattice data [16]. 

In Figs. |6| we show the results for A(q 2 ) and F(q 2 ) renormalized at three different points, 
IX = {4.3, 10,22} GeV with a(fL 2 ) = {0.21,0.16,0.13} respectively. On the right panel we plot 
the corresponding F(q 2 ) renormalized at the same points. Notice that the solutions obtained are in 
qualitative agreement with recent results from large-volume lattices [|TJ] where the both quantities, 
A(q ) wAF(q 2 ), reach finite (non- vanishing) values in the deep IR. 

The results for 1 + G(q 2 ) and L(q 2 ), renormalized at the same points, are presented in Fig. [7]. 
As we can see, the function 1 + G(q 2 ) is also IR finite exhibiting a plateau for values of q 2 < 0. lGeV 2 . 
In the UV region, we instead recover the perturbative behavior (2.14). On the other hand, L(q 2 ) 
shows a maximum in the intermediate momentum region, while, as expected, L(0) = 0. 
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— — — — — — | 0,8-| — — — — — — — 

1E-4 1E-3 0,01 0,1 1 10 100 1000 1E-4 1E-3 0,01 0,1 1 10 100 1000 

q 2 [GeV 2 ] q 2 [GeV 2 ] 

Figure 6: Left panel: Numerical solutions for the gluon propagator obtained from the SDE renormalized 
at three different points, n = {4.3, 10,22} GeV with a(jli 2 ) = {0.21,0.16,0.13}. Right panel: The ghost 
dressing function F(q 2 ) obtained from its corresponding SDE and renormalized at the same points. 



O 0,9 




q [GeV ] 




otdO = 0.21 and p = 4.3 GeV 
afti ! ) = 0.16 and (i= 10 GeV 
a(n') = 0.13 and n = 22 GeV 



0,000 

1000 1E-4 1E-3 0,01 0,1 



10 100 1000 10000 



q [GeV 2 ] 



Figure 7: Left panel: 1 +G(q 2 ) determined from Eq. (3.8), using the solutions for A(q 2 ) and D{q 2 ) pre- 
sented in the Fig. 6 at the same renormalization points. Right panel: The function L(q 2 ) obtained from 
Eq. (3.8). 



With all ingredients defined, the first thing one can check is that indeed Eq. ( [2.1 5[ ) gives rise to 
a RG-invariant combination. Using the latter definition, we can combine the different data sets for 
A(<7 ) and [1 + G(q 2 )] 2 at different renormalization points, to arrive at the curves shown on the left 
panel of Fig. |[ Indeed, we see that all curves, for different values of fL, merge one into the other 
proving that the combination d(q 2 ) is independent of the renormalization point chosen. 

From the dimensionful d(q 2 ) we must now extract a dimensionless factor, g 2 (q 2 ), correspond- 
ing to the running coupling (effective charge). Given that A(q 2 ) is IR finite, the physically mean- 
ingful procedure is to factor out from d{q 2 ) a massive propagator [q 2 + m 2 (q 2 )]~ l , 

zT^hs ' (4 - 4) 

q L +m z {q z ) 
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where for the mass we will assume "power-law" running JT7|], m 2 (q 2 ) =mf ) / (q 2 + ml). 

Thus, it follows from Eq. ([O]), that the effective charge a PT (q 2 ) = g 2 (q 2 )/4n is identified as 
being 

4KCC PT (q 2 ) = [q 2 + m 2 (q 2 )]d(q 2 ) . (4.5) 

Finally we compare numerically the two effective charges, a PT (^ 2 ) and ot g h(<7 ) on the right 
panel of Fig. ||. First, we determine a PT (q 2 ) obtained using ( p3| ), then we obtain a g h(q 2 ) with 
help of (3. 12) and the results for 1 + G(q 2 ) and L(q 2 ) shown in Fig. [7|. As we can clearly see, both 



couplings freeze at the same finite value, exhibiting a plateau for values of q 2 < 0.02 GeV 2 , while in 
the UV both show the expected perturbative behavior. They differ only slightly in the intermediate 
region where the values of L(q 2 ) are appreciable. 



5. Conclusions 

In this talk we have compared the definition of two QCD effective charges, a vl {q 2 ) and 
a g h(g 2 ), obtained within two vastly different frameworks: the PT-BFM on the one hand, and the 
ghost-gluon vertex (with the Taylor-kinematics) on the other. 

Despite their distinct field-theoretic origin, their dynamics involves the gluon propagator A(q 2 ) 
as a common ingredient and two different ingredients, which participate in a non-trivial identity. 
This identity, which is valid only in the Landau gauge, relates the ghost dressing function, F(q 2 ), 
with the two form-factors, G(q 2 ) and L(q 2 ). 

As consequence of the aforementioned identity, we have shown that the two effective charges 
are almost identical in the entire range of physical momenta. More specifically, they coincide ex- 
actly in the deep infrared, where they freeze at a common finite value, signaling the appearance of 
IR fixed point and a conformal window in QCD [Bt], in agreement with a variety of phenomenolog- 



ical studies [18]. 
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